Vortex Density of States and Absorption in Clean Layered Superconductors. 
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We study the spectrum of the states localized in the vortex cores in the mixed state of clean 
layered superconductors. S-wave coupling is assumed. It is found that in a large region of parameters 
adjacent to the superclean case a new universal (i.e. independent of the density of impurities) class 
of level statistics arises. It is the circular unitary random matrix ensemble. The density of states for 
such conditions is calculated. The absorption resulting from the Landau-Zener transitions between 
these levels is different from the classical result for an isotropic three-dimensional system. 
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Vortex cores in superconductors provide the low en- 
ergy states within the superconducting gap. These states 
are therefore important for the. absorption of low fre- 
quency electromagnetic field.ETu Of particular interest, 
are such states in layered s-wave superconductors,0Otil 
that mimic one of the features of high-T c superconduc- 
tors (HTSC) - their large anisotropy. As tunneling be- 
tween layers is slow, such systems can be considered quasi 
two-dimensional (2D). 



FIG. 1. DOS near the Fermi level in clean system. The 
result for class C of the,the random-matrix theory is shown 
by the dashed line. u'Q Eq. (^) is shown by the solid line. 
The figure also displays the result of the numerical so- 
lution of the system (H) (see below) averaged over 3000 
realizations of disorder (jl|), for 6 = 1, n^ 2 = 4 (four im- 
purities per vortex core), £/Af = 50. In the calculation 
we have used 100 basis CDM states. 

In the absence of disorder the spectrum of the Caroli 
- de Gennes - Matricon (CDM) states localized in a 2D 
vortex is equidistant: E® = —u)^{n— 1/2), loq ~ l/mf; 2 ~ 
A 2 /E-p, where n is an integer; £, m, A, E-p are the su- 
perconducting correlation length, mass of electron, su- 
perconducting gap, and Fermi energy respectively. Here 
and below we assume ft = 1, and measure energies from 
the Fermi level. In the presence of a small number of 
the short-range impurities the spectrum retains some of 
its periodicity. In particular in the superclean regime 
(1/r <C u>o) such impurities displace the levels with 
even and odd n in the opposite directions by an equal 



amount.B This implies that the spectrum is 2loq periodic. 
This result holds if the impurities are strong, i.e. their 

1/2 

Born parameter 9 » (A/E F ) 1 . 

The levels are shifted only if an impurity is present in- 
side the vortex core. Because in the superclean regime 
the density of impurities is small such vortices are rare. 
Hence, the shape of DOS in the superclean regime de- 
pends on the density of impurities and is not universal. 

As the amount of disorder is increased the clean regime 
(ujq <C 1/t <C A) is approached. In this regime the 
vortex cores can be viewed simpJistically as small gran- 
ules of disordered normal metal.u Therefore one can_£x=. 
pect the energy level statistics of Wigner-Dyson typeB'EHl 
More precisely, due to the superconductor surrounding 
the "granule" , one can expect level statistics of class C 
in the nomenclature of Ref [?]. This is indeed the case 
as shown in Ref. |^, for layered superconductors with a 
white noise disorder potential, i.e. for vanishingly small 
Born parameter 9. The corresponding density of states 
(DOS) is shown in Fig. |l| by a dashed line. The suppres- 
sion of DOS on the Fermi level is due to the particle-hole 
symmetry (E — > —E) pertinent to the quasiparticles. 

In this work we study the clean superconductors with 
strong short-range impurities: 



Vi m p(r) 



Vi5 (r - Oi) , 
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where Vi and are the strength and the position of 
the z-th impurity respectively, while S(r) is the func- 
tion with spatial extent comparable to Ap. We assume 
that the impurity scattering phase 9 — miy 2 ) 1 ! 2 is of 
the order of one. We address the question of whether 
the 2u>o periodicity of the spectrum, specific for the 
superclean regime, persists in the clean regime. We 
show that in the large subregion of the latter defined 
by 9 2 AyJ A/Ep » 1/t 3> wot the spectrum is indeed 
2ujq periodic. Thus it is effectively defined by positions 
of the two nearest levels, with all the others multiplied 
periodically up and down. The DOS in our case is given 
by 



P(E) 
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It is shown in Fig. |l| by solid line. Due to the particle- 
hole symmetry of the spectrum and the periodicity, the 
average DOS also defines the level spacing distribution 
p(s) through the relation p(s) = p(s/2)/2, where s = 2E 
is the level spacing. The distribution (0) is then consist 
tent with the unitary circular random matrix ensemblcEl 
for the two independent levels. This is in contrast to the 
findings of Ref. ^ where white noise disorder is consid- 
ered (9 — > 0). We conclude therefore that in the above 
limit a new universal (i. e. independent of the density of 
impurities) class of level statistics arises. 

We calculate also the dissipative componant-pf conduc- 
tivity due to the Landau-Zener transitions :□□ 



<r xx = 32a n H c2 lnO 1 /B(uj t) 



(3) 



Because in the clean case the parameter luqt <C 1 this re- 
sult obtained for the layered superconductocs is different 
from the isotropic three-dimensional caseffij 



Here (oj,^) is the position of the i-th impurity in the 
polar coordinates. In Eq. (||) we ignore the influence of 
impurities on the order parameter. The summation index 
to in Eq. (|J) runs from — N to N w A/ojq, where N is a 
very large number. We will assume further N to be infi- 
nite and perform a Fourier transform from integer index 
to to the real variable <j>, usine 
After this transformation Eq. (\ 
form 



-i(pm 



I) acquires the following 



- ioj ^ + j~" L(<p, <f>'M<f>'W = (e~y) *M> 

(10) 

where the kernel L(<p,(f)') describes scattering by impuri- 
ties: 

W) =E^ A ^[ 1 - ex p(^-^)] 

i 

cxp (— ikpdi [sin((/> + 4>i) — sin(</)' + 4>i)\) 



2tt 



0~xx ~ <J n H c 2/ B. 



(4) 



The two-component wavefunctions of excitations ip 
(ipi = Uw—M>2 — v) satisfy Bogolubov - de Genncs 
equations Jla 



[<Jz (H + V imp ) + a x ReA(r) + a y lmA(r)} i/j = Exp. (5) 

Here Hq — p 2 /2m — Ep, {a x , a y , a z } are the Pauli ma- 
trices, and E is the excitation energy. In the absence 
of impurities the wavefunction of the n-th CDM state is 
given by 



. e i(n-l)0 j n _ % (fcpr) 



(6) 



Here J n (x) is the Bessel function and K(r) — Ar/wp = 
r/^n. Here and below we assume that the absolute value 
of the order parameter is constant within the vortex 
core at|-small temperatures, adopting the Kramer-Pesch 
ansatz.tHl If impurities are present the excitation wave- 
function can be found as a mixture of different CDM 
states with coefficients C m : 



(7) 



C m satisfy the general perturbation theory equation 

(E°-E)C n +J2KmCm=0, (8) 



where A l nm is the matrix element produced by one short- 
range impurity between two CDM states n and m. The 
expression for this matrix element follows from Eqs. (Q), 
(§, and (1): 

Km = V ie - 2K ^ +i ^ m -^ [J n (fc F a,) J m (k F a t ) (Q) 

-J„_l (kpOi) J m -l (kpOi)] Af£- 



FIG. 2. Motion of a quasiparticle within the vortex 
core. The center of the core is shown by the dark region. 
The solid line shows a piece of the trajectory of the quasi- 
particle. The arrow indicates the direction of momentum 
of the quasiparticle, which is unchanged by the Andreev 
reflection. Impurity scattering however changes it to the 
opposite. The impurities are shown by the black spots. 
When colliding with an impurity the quasiparticle scat- 
ters into the state, which is almost diametrically opposite 
(dashed line). 

The physical picture behind this equation is as fol- 
lows: the quasiparticle flies back and forth almost in a 
straight Jine, Andreev reflecting from the walls of the vor- 
tex coretj (see Fig. |l|). This motion is performed with 
the characteristic frequency i>f/£ ~ A, which is relatively 
large. In addition to this the quasiparticle performs a 
slow precession with the frequency wo « A, described 
by the first term in Eq. (|10|) . 

The second term in Eq. (|10|) describes scattering of the 
quasiparticles by impurities. n Let us first consider a core 
with only one impurity in it.H It is typically situated at a 
distance a, ~ £ from the center of the core. In this case 
the parameter k^cn is large. The scattering matrix (111]) 
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is then localized at the angles fa ± tt/2, with the charac- 
teristic spread 5 fa ~ 1/y/k-pai <C 1. This means that the 
impurity scatters the quasiparticle from d> = fa + tt/2 to 
</), — 7r/2 and vice versa (see Fig. ||). Eq. pQ ) can then be 
integrated near these points to render the boundary con- 
ditions for the wavefunction at them. To present these 
boundary conditions it is convenient to introduce instead 
of one function *5/(<fi) the spinor 



' (0 + 7T) 



which is naturally defined for < 
boundary conditions read: 



(12) 



< 7r. Then the 



* (fa + n/2 + v) = (fa + tt/2 - v) , (13) 

where v 3> Solving Eqs. (|l0| ) and ( |IT|) near impurity 
we obtain the transfer matrix at the impurity 



Mi = 



1 



1 - |jf -2J 
2J* 1-|J| : 



(14) 



where J = V;e 2ifepa4 - 2i ^ ai )/27ru;oai£- Noticing that be- 
tween collisions with impurities according to Eq. (|l^ 

i& oc exp [i (E/ujq — 1/2) fa , and the boundary condi- 
tions for spinor at points = and tt are ^(-n - — 0) = 
cr x $(+0), we obtain the equation for the excitation en 
ergy levels in the presence of one impurity 



cos(itE/uj ) = -iTi(a x Mi) /2. 



(15) 



This equation coincides with the one derived earlier in 
Ref. |. 

This approach can be generalized to include many im- 
purities. In this case the transfer matrices corresponding 
to different impurities should be multiplied to provide the 
transfer matrix from (f> = to tt. The resulting equation 
for the energy levels inside the vortex core is analogous 
to (p|) 



cos (nE/wo) = -iTr ( a x J] M t /2. 



(16) 



This equation is the central result of this paper. It is 
valid in the limit Ap — > 0. More precisely NiSfa ~ 
Ni(A F /£) 1/2 < !■ Here Ni ~ n^ 2 is the number of 
impurities in the vortex core and n, is the density of im- 
purities. 

From Eq. (|l^) it follows that the spectrum of the sys- 
tem is periodic with the period 2loq. It is also symmetric 
with respect to the replacement E — >• —E (particle-hole 
symmetry). Let us study the properties of this spec- 
trum in the limit Ni 3> 1. The unitary matrix Mi is 
a quaternion. Therefore it can be written in the form 
Mi = noicro + iriicr, where uo is the unit 2 by 2 matrix, 
and rii = (noi,rii) is a four-dimensional (4D) unit vec- 
tor. The product of the quaternions is also a quaternion. 
Therefore 



M = Y[ Mi - riQ<TQ + incr, 

i 

cos (ttE/uq) = n x , 



(17) 



where the unit 4- vector n parameterize the matrix M. 
Every impurity produces a rotation of this vector. After 
the action of many impurities the probability to find it at 
a particular point on 4D sphere n 2 = 1 spreads uniformly 
over the 3D surface of this sphere. Therefore averaging 
over the position of impurities can be rep laced by aver- 
aging over the 4D unit sphere. Eq. (|l6|) results in the 
following expression for the average DOS 



P(E) 



d 4 n5 (n 2 — l) 8 [E — cj arccos (n x ) /tt] (18) 



This results in Eq. (S) above. 

We now turn to the calculation of vortex viscosity in 
this model. The vortex motion produces transitions be- 
tween excited states within the vortex core. These transi- 
tions result in the energy dissipation equal to loqW, where 
W is the transition rate between nearest levels.li-3 Viscos- 
ity r) can be related to the dissipation rate per vortex 
core Q through: 



r/ = Q/v 2 = Lu W/v 2 



(19) 



Here v is the velocity of vortex motion. Our task will be 
to determine W. We assume that the velocity is small, 
so that the J^naitions between levels are of the Landau- 
Zener type.H&oEj Levels therefore have to come very 
close to each other. Due to the particle-hole symmetry 
of the spectrum this implies that energy of the states near 
Fermi level should be very close to zero. From Eq. ( |l7j ) 
we see that in this case the component n x of the 4- vector 
identifying the position of matrix M on the 4D sphere 
should be close to 1. Let us imagine now that the vortex 
is moving relative to the impurities. Vector n will then 
wander over the surface of the 4D sphere occasionally ap- 
proaching the pole n x = 1, where the Zener transitions 
are possible. Due to rotational symmetry we can assume 
that it approaches the pole from the direction parallel to 
no. Therefore no changes sign near the pole and can be 
written as hq — v^t, where i>4 = n is the 4-velocity of 
the point on the sphere. Expanding the expression for 
energy in Eq. (]lj) we obtain the distance between levels 
when they are close to each other as a function of time: 



s(t) = 2uj a ^Jn 2 y + n 2 z + vlt 2 /TT. 



(20) 



This allows us to use the Landau-Zener formulaE! to cal- 
culate the probability of transitions: 



W = dn y dn z V4 exp 



^. (21) 

UJQ 
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Here dn y dn z V4 is the differential frequency of attempts. 
We conclude that the average rate of Landau-Zener tran- 
sition is proportional to the average square of the velocity 
of vector n. Our next goal will be to evaluate this aver- 
age. 

Vector n moves on the 4D sphere because the matri- 
ces Mi in Eq. ([l7]) change. The matrices change due to 
the vortex motion relative to the impurities, according 
to Eq. (14). Noticing that the average velocity V4 is the 
same all over the sphere because the density of vector 
11 is uniform there, and that all 4 directions of V4 are 
equivalent we write 



(22) 



The average is taken over the impurity potential. The 
latter equation emphasizes that the average of the cross 
term containing the products of two impurities is zero. 
To evaluate the derivative with respect to the coordinate 
of one impurity dng jdxi we notice from Eq. ( [l7|) that 



n = Tr (MiR) /2. 



(23) 



Here R is a random 2 by 2 unitary matrix obtained from 
M by cyclic permutations under the trace. Considering 
then a small displacement of the impurity bxi we obtain 



Im(J Ru} xj 

dn = 4fc F _. I -7-10 oxi 

1 + \J\ 2 at 



(24) 



We then average (Shq/ 'Sxi) 2 over the random matrix R 
instead of impurities: 



,,2 _ 



2irrdr 



iv 2 k 2 \J(r)\ 



1 2-KV 2 kl, 1 



(i+m 2 r 



LUQT 



In-. (25) 



Using this result, Eqs. (|19[), (|2l|), and the connection be- 
tween viscosity and conductivity a xx = \e\cq/wB, where 
e is the charge of electron, we obtain Eq. (||) above. 

In our consideration we ignored the energy relaxation 
processes, assuming the inelastic collision time r 6 to be 
large. More precisely, we have assumed that v 3> v e ~ 
Af-v^two/te. In the opposite case v <C v £ , Eq. (||) 
should be multiplied by a small factor v/v £ . At large 
velocities our approximation fails when the transitions 
to the next nearest neighbor levels become possible, i.e. 



v > Ap-i 



Our approach therefore has an interval of 



validity if wot £ ^> 1. We have also ignored the overheat- 
ing of the core, assuming r e v 2 <§; tA|A 2 , and pinning, 
assuming that either the observations are, made at a fre- 
quency exceeding the pinning frequency^ or the vortex 
lattice is rigid. 

The well-known examples of clean layered supietcon- 
ductors ace HTSC. For them both the vortex DOSE£l and 
viscosityiiH are measured. Due to the nodes in the super- 
conducting gap and the absence of the Fermi spectrum in 
the directions distant from the nodes, even with no impu- 
rities the spectrum of vortex states in HTSC is not clear. 
The applicability of our results to this case is not obvi- 
ous and warrants further study. We believe however that 



our conclusion about the existence of the circular random 
matrix ensemble is likely to be true for the d-wave super- 
conductors. This is justified by the recent study of clean, 
d-wave superconductors, revealing the 2ujq periodicity.E£l 
We believe therefore that DOS given by Eq. (0) may be 
observed, should our regime be accessible experimentally. 

In conclusion, we have studied the DOS and absorption 
by the vortex core states in clean layered s-wave super- 
conductors, when 8 2 Ay/ A/Ef > 1/t > uo a . We predict 
the existence of the circular random matrix ensemble in 
this regime. We have found that both DOS and the dis- 
sipative component of conductivity are different from the 
isotropic 3D case. This difference is due to the fact that 
in 2D the quasiparticles can scatter coherently on the 
same impurity profile many times. 
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